A free boundary problem arising from DCIS mathematical model

Ductal carcinoma in situ – a special cancer – is confined within the breast ductal only. We derive the mathematical ductal carcinoma in situ model in a form of a nonlinear parabolic equation with initial, boundary, and free boundary conditions. Existence, uniqueness, and stability of problem are proved. Algorithm and illustrative examples are included to demonstrate the validity and applicability of the technique in this paper. Copyright © 2016 John Wiley & Sons, Ltd.     

Study on convergence and error of a numerical method for solving systems of nonlinear Fredholm–Volterra integral equations of Hammerstein type

The aim of this paper is to present an efficient analytical and numerical procedure for solving systems of nonlinear Fredholm–Volterra integral equations of the Hammerstein type with the aid of fixed point techniques and the usual Schauder basis in an adequate Banach space.     

A Buzano type inequality for two Hermitian forms and applications

A Buzano type inequality for two nonnegative Hermitian forms is obtained. Applications to inequalities for norm and numerical radius of bounded linear operators in complex Hilbert spaces are given.     

Finite element analysis and approximation of Burgers’‐Fisher equation

In this article, the authors present finite element analysis and approximation of Burgers’‐Fisher equation. Existence and uniqueness of weak solution is proved by Galerkin's finite element method for non‐smooth initial data. Next, a priori error estimates of semi‐discrete solution in norm, are derived and the convergence of semi‐discrete solution is established. Then, fully discretization of the problem is done with the help of Euler's backward method. The nonlinearity is removed by lagging it to previous known level. The scheme is found to be convergent. Positivity of fully discrete solution is discussed, and bounds on time step are discovered for which the solution preserves its positivity. Finally, numerical experiments are performed on some examples to demonstrate the effectiveness of the scheme. The proposed scheme found to be fast, easy and accurate.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1652–1677, 2017     

An inverse backward problem for degenerate parabolic equations

This work studies the inverse problem of reconstructing an initial value function in the degenerate parabolic equation using the final measurement data. Problems of this type have important applications in the field of financial engineering. Being different from other inverse backward parabolic problems, the mathematical model in our article may be allowed to degenerate at some part of boundaries, which may lead to the corresponding boundary conditions missing. The conditional stability of the solution is obtained using the logarithmic convexity method. A finite difference scheme is constructed to solve the direct problem and the corresponding stability and convergence are proved. The Landweber iteration algorithm is applied to the inverse problem and some typical numerical experiments are also performed in the paper. The numerical results show that the proposed method is stable and the unknown initial value is recovered very well.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1900–1923, 2017     

Convergence analysis for second‐order accurate schemes for the periodic nonlocal Allen‐Cahn and Cahn‐Hilliard equations

In this paper, we provide a detailed convergence analysis for fully discrete second‐order (in both time and space) numerical schemes for nonlocal Allen‐Cahn and nonlocal Cahn‐Hilliard equations. The unconditional unique solvability and energy stability ensures ? ⁴ stability. The convergence analysis for the nonlocal Allen‐Cahn equation follows the standard procedure of consistency and stability estimate for the numerical error function. For the nonlocal Cahn‐Hilliard equation, because of the complicated form of the nonlinear term, a careful expansion of its discrete gradient is undertaken, and an H ^?1 inner‐product estimate of this nonlinear numerical error is derived to establish convergence. In addition, an a priori bound of the numerical solution at the discrete level is needed in the error estimate. Such a bound can be obtained by performing a higher order consistency analysis by using asymptotic expansions for the numerical solution. Following the technique originally proposed by Strang (eg, 1964), instead of the standard comparison between the exact and numerical solutions, an error estimate between the numerical solution and the constructed approximate solution yields an O (s ³+h ⁴) convergence in norm, in which s and h denote the time step and spatial mesh sizes, respectively. This in turn leads to the necessary bound under a standard constraint s ≤C h . Here, we also prove convergence of the scheme in the maximum norm under the same constraint.     

Quenching phenomenon of a time‐fractional diffusion equation with singular source term

In this paper, a time‐fractional diffusion equation with singular source term is considered. The Caputo fractional derivative with order 0<α ?1 is applied to the temporal variable. Under specific initial and boundary conditions, we find that the time‐fractional diffusion equation presents quenching solution that is not globally well‐defined as time goes to infinity. The quenching time is estimated by using the eigenfunction of linear fractional diffusion equation. Moreover, by implementing a finite difference scheme, we give some numerical simulations to demonstrate the theoretical analysis. Copyright © 2017 John Wiley & Sons, Ltd.     

一种新的求非线性方程组的数值延拓法

针对迭代过程中的Jacobi奇异问题,本文提出了一种新的数值延拓法.通过构造双参数同伦算子,采用可控条件和适当选取参数的方式克服Jacobi奇异性,并分析了方法的收敛性.最后,通过数值实验对比,验证了方法的可行性和优越性.特别是具有可调控越过Jacobi奇异(点、线、面)的优势,从而也在某种程度上解决了数值延拓法严重依赖于初值的问题.     

分布式通信响应优化问题及其内点法求解

为了优化移动互联网的分布式通信系统的响应速度,建立分布式通信系统响应速度最优化问题的数学模型,并设计和改进求解该最优化问题的内点法.针对该最优化问题发展一套高效率预条件方法来帮助求解内点法,不但改善计算方法的数值稳定性,而且提高算法的计算效率.通过数值实验验证该预条件对算法稳定性和效率的提高.     

Variational Bayes With Intractable Likelihood

Variational Bayes (VB) is rapidly becoming a popular tool for Bayesian inference in statistical modeling. However, the existing VB algorithms are restricted to cases where the likelihood is tractable, which precludes their use in many interesting situations such as in state--space models and in approximate Bayesian computation (ABC), where application of VB methods was previously impossible. This article extends the scope of application of VB to cases where the likelihood is intractable, but can be estimated unbiasedly. The proposed VB method therefore makes it possible to carry out Bayesian inference in many statistical applications, including state--space models and ABC. The method is generic in the sense that it can be applied to almost all statistical models without requiring too much model-based derivation, which is a drawback of many existing VB algorithms. We also show how the proposed method can be used to obtain highly accurate VB approximations of marginal posterior distributions. Supplementary material for this article is available online.